Theory of Aging in Structural Glasses 
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The random first order transition theory of the dynamics of supercooled liquids is extended to 
treat aging phenomena in nonequilibrium structural glasses. A reformulation of the idea of "entropic 
droplets" in terms of libraries of local energy landscapes is introduced which treats in a uniform way 
the supercooled liquid (reproducing earlier results) and glassy regimes. The resulting microscopic 
theory of aging makes contact with the Nayaranaswamy-Moynihan-Tool nonlinear relaxation for- 
malism and the Hodge-Scherer extrapolation of the Adam-Gibbs formula, but deviations from both 
approaches are predicted and shown to be consistent with experiment. The nonlinearity of glassy 
relaxation is shown to quantitatively correlate with liquid fragility. The residual nonArrhenius 
temperature dependence of relaxation observed in quenched glasses is explained. The broaden- 
ing of relaxation spectra in the nonequilibrium glass with decreasing temperature is quantitatively 
predicted. The theory leads to the prediction of spatially fluctuating Active temperatures in the long- 
aged glassy state, which have non-Gaussian statistics. This can give rise to "ultra-slow" relaxations 
in systems after deep quenches. 

PACS numbers: 



I. INTRODUCTION 

The energy landscape metaphor has turned theorists 
towards viewing the global geometry of the phase space 
of complex systems. When these systems are mesoscopic 
in size, for example, small proteins or gas phase clus- 
ters [2j, a more or less complete mathematical formula- 
tion of the idea, capable of treating kinetics and ther- 
modynamics, can be made using statistical tools to char- 
acterize minima and saddle points of the entire system. 
Yet for macroscopic systems, most transitions rearrange 
particles only locally. This essential aspect of the dy- 
namics is brought home forcefully by noting that a liter 
of liquid water will move from one energy minimum to 
another in 1CP 39 sec. Such a short time scale cannot 
be directly relevant to any laboratory measurement on 
this system! The necessity for using a local description 
is clearly recognized in the modern theory of supercooled 
liquids and glasses which is based on the statistical me- 
chanics of random first order transitions 0-0 • In the 
deeply supercooled liquid regime this theory explains not 
only the phenomenological features of the dynamics 
but quantitatively predicts, on a microscopic basis, the 
size of cooperative lengths, the precise non-Arrhenius be- 
havior of typical relaxation times and the non- 
exponentiality of relaxation 0. A quantized version of 
the theory explains the low temperature thermodynamics 
of amorphous substances, usually interpreted in terms of 
two level systems and the more energetic Boson peak 
excitations The crucial manifestation of the local- 

ity concept in this theory (which has many mean field, 
global, aspects) is the notion that "entropic droplets" 
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drive the large-scale activated notions in glass forming 
liquids |J- LJ li| and give both liquids and glasses an in- 
trinsic "dynamical mosaic" structure. 

The main purpose of this paper is to describe the pre- 
dictions of the RFOT theory for the behavior of glasses 
that have fallen out of equilibrium because of being 
rapidly cooled from the melt. Dynamics does not cease in 
rapidly quenched liquids that become glasses; rather mo- 
tions persist but are frustratingly slow and hard to study 
experimentally. This slow, far-from-equilibrium dynam- 
ics, called "aging", is not only of fundamental interest 
for statistical mechanics ^lj but also is of great practical 
interest since so many amorphous materials are used in 
everyday life for times exceeding their preparation time 
|12|. Sometimes even small changes of properties upon 
aging are crucial to engineering performance. 

Like the quantum theory of the low temperature prop- 
erties of glasses, the description of the nonequilibrium ag- 
ing regime requires the explicit construction and study of 
the energy landscape of local regions of the glass. This lo- 
cal energy landscape description turns out to be a micro- 
canonical ensemble reformulation of the entropic droplet 
concept that was formulated originally in the canonical 
ensemble. We have found this microcanonical description 
to be more vivid and somewhat easier to communicate 
than the original canonical formulation. After introduc- 
ing this formulation we will show that it indeed repro- 
duces the results for the dynamics of equilibrated super- 
cooled liquids already obtained using the random first or- 
der transition theory. More important, this local energy 
landscape theory makes several very striking predictions 
about the aging regime and how the aging dynamics in 
the nonequilibrium glass is related to the kinetics and 
thermodynamics of the equilibrated liquid. The predic- 
tions are consonant with all experiments known to us 
and make contact with, but are formally distinct from, 
the phenomenological approaches commonly used to de- 
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scribe aging in structural glasses [H M M H 113 The 
present theory, for example, makes a specific prediction 
of the so-called nonlinearity parameter, Xnmt, in the 
Narayanaswamy-Moynihan-Tool, formalism [131 fl4L fl5| . 
In addition, the degree of nonexponentiality of relaxation 
characterized by a (3 exponent is predicted. turns out to 
be neither precisely fixed at T g nor does it precisely scale 
like the (3 for an equilibrated liquid with T, another as- 
sumption often made. The predicted variation of (3 with 
quench temperature is however modest until very low 
temperatures are reached. Likewise the Adam-Gibbs ex- 
pression [l8j has been extended phenomenologically no- 
tably by Hodge and Scherer into the aging regime below 
T g by assuming the configurational entropy to be frozen 
at T g . While above T g , an Adam-Gibbs like form of the 
temperature dependence of relaxation time is found from 
the present theory, below T g , the local energy landscape 
approach gives a different expression for the relaxation 
rate in an aging glass. Quantitatively this expression 
gives relaxation rates close to the Hodge-Scherer- Adam- 
Gibbs (HSAG) latter extrapolation but the local energy 
landscape theory predicts a deviation from that formula. 
This deviation would be interpreted within the HSAG 
framework as a quench temperature dependence of the 
apparent configurational entropy. Such a deviation has 
been found by Alegria et al. in their careful work on 
aging of polymers |19|. 

Our most explicit results are obtained for the idealized 
situation of the aging initially found after a cooling his- 
tory with a single rapid quench of modest depth. We will 
also qualitatively discuss modifications of the simple the- 
ory expected for very deep quenches. We also discuss the 
behavior of systems that have significantly relaxed in the 
quenched state. While, our analysis suggests that, to a 
first approximation, introducing a single Active tempera- 
ture should serve well to describe many quench histories, 
the local energy landscape theory shows that using a sin- 
gle fictive temperature is not exact. We suggest possible 
modifications of the usual aging kinematics based on the 
present theory. 

The organization of the paper is as follows: we first 
describe how the local energy landscape view of entropic 
droplets can be visualized and show how a library of lo- 
cal energy landscapes can be constructed. We next show 
how to derive in this framework the (previously obtained) 
Vogel-Fulcher behavior above T g . We then present the 
results for typical relaxation rate for rapidly quenched, 
aging glasses and compare these predictions with exper- 
imental results. We then discuss the predictions for the 
stretching exponent both above and below T g . Following 
this, the generalizations needed for very deeply quenched 
glasses and glasses that have undergone significant relax- 
ation in the aging state are discussed. Finally we sum- 
marize our theory and suggest some further experimental 
tests of it. 



II. THE LOCAL ENERGY LANDSCAPE 
CONSTRUCTION 

The energy landscape language is usually applied to a 
small system (protein or cluster). It is also used for a 
thermodynamically large system described by mean field 
theory [20L l2*lj . In the first situation, the barriers be- 
tween states are, of course, finite because the system is 
finite, but barriers may be formally infinite for the mean 
field system. The "states" for a cluster or a protein are 
often taken to be the basins surrounding minima of the 
potential energy |22( . These are well defined and no tran- 
sitions can occur, classically, between them at absolute 
zero. The "states" of a mean field system are tied to min- 
ima of free energy and again, owing to the O(N) barriers, 
no transition can occur between them (even at finite T). 

In a supercooled liquid the observed plateau in the 
time dependent neutron scattering correlation function 
|23| shows that most molecules spend a great deal of time 
vibrating about a given location. At least deep in the su- 
percooled regime (where the plateau is well developed) 
we can, therefore, conceptually imagine constructing an 
average location for any particle about which it vibrates 
(for times less than the plateau). The three dimensional 
structure based on these average locations will be quite 
close to a potential energy minimum found by remov- 
ing kinetic energy from the system with infinite speed 
( "steepest descent" to an "inherent structure" ) . A mean 
field theory of the glass transition can be obtained by 
constructing a free energy functional dependent on the 
(plateau time averaged) density and then noting that 
this functional has minima for densi ty p atterns localized 
around such inherent structures 0> EH ■ I n this way 
a vibrational component of the free energy can be de- 
fined and can be added to the average energy of such a 
state to give a free energy which determines the thermal 
probability of being in this state through the Boltzmann 
law. 

Free energy functionals only have such aperiodic min- 
ima below a temperature Ta- This temperature has been 
shown (within a simple approximation) to be equivalent 
to the mode coupling dynamical transition temperature 
|3|. The transition has the character of a spinodal for 
a first order phase transition [2J. Above Ta even simple 
vibrational motions will take the system from one mini- 
mum to another, but below Ta, the plateau in the time 
dependent structure factor indicates the persistence of 
such states and therefore their relevance to dynamics. In 
fact, as a first approximation, the persistence of these 
states allows them to act as a "basis set" for describing 
supercooled liquid dynamics. Thus we say a free energy 
landscape emerges at Ta- Ta has been evaluated both for 
simple liquids |24i |25| and for models of network form- 
ing liquids with repulsive force p?) from first principles. 
The physical meaning of the temperature Ta is simple. 
When supecooled below its melting temperature T m , the 
liquid is, of course, in a metastable state as a whole, and 
therefore must be metastable locally. Ta can be conve- 



niently interpreted, for example, as a temperature below 
which every two molecules will spend a particular num- 
ber (say, 300) of vibrational periods together between the 
first encounter and the final departure. 

As part of our programme to describe the microscop- 
ies of liquid relaxations, we will show that transitions 
between the metastable configurations of the liquid as 
a whole consists of transitions between local metastable 
configurations. We will employ a microcanonical proce- 
dure in which only local regions are considered and will 
find that beyond a certain (relatively small) size N* , the 
thermal and relaxational properties of the liquid do not 
depend on the size of the sample. Alternatively said, all 
liquid properties of interest can be deduced by focusing 
on a liquid region of size N* and completely disregarding 
what the rest of the liquid is doing. This is the essence 
of the locality of the liquid free energy landscape. The 
microcanonical procedure is a necessary step in estab- 
lishing such locality; it is however rarely used, so let us 
first train our intuition on the very familiar example of 
a harmonic lattice. Imagine being inside an extremely 
large, cubic, and purely harmonic lattice and being given 
the ability to do arbitrary thermodynamic measurement 
locally. You are further assigned to explain, within a for- 
mal model, those thermal measurements, but are allowed 
to visually inspect and mechanically test the bonds with 
an arbitrarily large, but finite region, limited by how fast 
you can perform the inspection. Upon checking that all 
individual bonds are truly harmonic and the lattice is 
cubic within a certain region (of size L), you write down 
a simple hamiltonian but are left with the issue of the 
conditions at the boundary - strictly speaking, the as- 
signement is undoable! Being an optimist, you say: let 
me assume for today that the rest of the sample is totally 
rigid (fixed boundary conditions) or does not exist at all 
(open bondary condition) and diagonalize the resulting 
Hamiltonian; tomorrow I can expand my horizon, repeat 
the procedure and see what happens. The first day proves 
frustrating though, because the Hamiltonian has an en- 
ergy gap (proportional to tt/L), while the thermal mea- 
surements clearly show very low frequency excitations 
are present. Fortunately, since (unbeknownst to you) 
the lattice was harmonic and periodic, your consecutive 
inspections and diagonalizations will yield a smaller and 
a smaller gap. Clearly, the landscape of a regular har- 
monic lattice is non-local and the microcanonical proce- 
dure offers a definitive test of non-locality. Consider yet 
another elementary notion: Suppose you have a prod- 
uct of two Hilbert spaces A ® B (possibly interacting 
via Va,b) with an energy £i A ,j B assigned to configura- 
tions ia and js of the subspaces A and B. Consider the 
partition function Yli A JB e - ^^ ^ = J2i A e - ^, where 
e -0<t>i A = e~P ei A-3B . If A and B are independent: 
Va.b = => £i A ,j B — e iA Jrt jB > ~ one ma y further simplify 
4>i A = ti A ,j B — TSb- In either case, tpi A 's are free energies 
of the degree of freedom B, but can still be (with advan- 
tage) regarded as an energy, as far as the A degree of free- 
dom is concerned. For example, one is allowed to build a 
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FIG. 1: In the upper panel on the left a global configuration 
is shown, chosen out of a global energy landscape. A region 
of N = 5 particles in this configuration is rearranged in the 
center illustration. The original particle positions are indi- 
cated with dashed lines. A larger rearranged region involving 
N = 7 particles is connected dynamically to these states and 
is shown on the right. In the lower panel, the left most fig- 
ure shows the huge density of states that is possible initially. 
The density of states found in the local library originating 
from a given initial state with 5 particles being allowed to 
move locally is shown in the second diagram. These ener- 
gies are generally higher than the original state owing to the 
mismatch at the borders. The larger density of states where 
7 particles are allowed to move is shown in the right most 
part of this panel. As the library grows in size, the states 
as a whole are still found at higher energies but the width of 
the distribution grows. Eventually with growing N, a state 
within thermal reach of the initial state will be found. At this 
value of N* we expect a region to be able to equilibrate. 



microcanonical construction using <\>i A 's as the label. In 
the following, we will use a 0^-like quantity to describe 
all the degree of freedom in the supercooled (or quenched) 
liquid in excess of the lowest energy crystalline structure 
corresponding to this substance (at this temperature). 
(For a polymer that does not form a crystal our results 
will still apply qualitatively.) While the "integrated out" 
degrees of freedom B are clearly related to the vibrations 
of the corresponding crystal structure, they are easy to 
conceptualize only in the mean field limit, where they are 
indeed harmonic vibrations (see shortly). Otherwise, the 
"£?" motions are strongly anharmonic. 

We are now ready construct a library of free energy 
minima for a very large sample of a liquid (or glass). 
This spectrum is shown schematically in the first column 
of Figure n as a function of the free energy of a state i, 
cj) l i lb which we take as the sum of an energy ej and an 
cntropic contribution from vibrations within the basin 

Because the system is very large this global library of 
states has an exceedingly dense spectrum, whose density 
dramatically increases with free energy. Let us imag- 
ine the system is presently in one of these basins with 
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(free) energy (jyf' , We can now construct a set of local 
libraries of states. To do this, imagine mentally cutting 
out a region around a location R containing N particles, 
where N is much less than the total number of particles 
in the sample. Call this region D(R,N). Next, freeze 
the molecules outside this region but allow molecules 
within D to move. (Unlike in the cubic crystal analogy, 
a supecooled structure has built-in stress, so using open 
boundary condition is not useful.) With the frozen en- 
vironment one could (by quenching the potential energy 
or a free energy functional) find a new set of inherent 
structures that involve only reconfiguring the particles in 
D. Actually such a set of structures (apart from strains 
near the surface of D) will locally resemble a subset of 
the original global library. If N is small, however, we ex- 
pect the subset of the states so sampled to be selected on 
average from a higher energy part of the global spectrum 
than where the original <p l i lb was found. Essentially this 
is because the density of states rapidly increases with 
</> and there is less freedom to readjust particles in the 
constrained minimization problem because of the fixed 
environment. In absolute terms, we expect the discrep- 
ancy in energy of the mean energy of the states in any 
local library from the starting state to increase with N 
since it is the particles at the borders that must be most 
strained. 

When N is small the local library is sparse but it grows 
denser and spreads out over a larger range of free energy 
as the size of the region increases. The set of such li- 
braries centered around R is shown also in Figure ^ for 
N = 1 up to a large value of N. The number of states 
in the library at free energy <pi determines the configu- 
rational entropy Q c (<j) l i lb ) — e Sc ^* where fee is the 

Boltzmann constant. The competition between the av- 
erage energy growth and the spreading of the range of 
energies with size of the region means that there is a 
characteristic size N* where a state will finally be found 
within thermal energy of the starting state. N* turns out 
to be the size of a dynamically correlated region in the 
liquid or glass. Since the new configuration is within fc^T 
of the starting state, a transition of such a region from 
its original configuration to the new state can occur with 
reasonable probability. The region after reconfiguration 
will be characterized by a temperature T. Elementary 
transitions in the liquid must leave most of the molecules 
near their old locations - this is the dynamical essence 
of locality. Therefore the local landscape libraries are lo- 
cally connected to each other: in order to re-arrange a 
large region, smaller regions located in the same place 
must first be rearranged. The specific law connecting 
states in neighboring libraries must obey detailed bal- 
ance (at the vibrational temperature T) but otherwise 
will depend specifically on molecular details for each sys- 
tem. The additional activation energy, for a downhill 
move in such a library should only be a few times the 
energy needed to heat a particle to Ta, the dynamical 
transition temperature. 

Irrespective of the detailed motions allowed, the local- 



ity of the dynamics, however, guarantees that the tran- 
sitions from the initial state to one of the thermally al- 
lowed states in the library for -D(R, N*) will nevertheless 
be slow: owing to the initial rise of the average energy 
there will be a bottleneck in the probability flux at an 
intermediate value of N, namely N$. The activation free 
energy for reaching this bottleneck state determines the 
escape rate from the initial configuration by motions in 
the vicinity of R. 

While the local libraries can be constructed explicitly 
given sufficient computational resources, the RFOT the- 
ory suggests a useful set of approximations to the statis- 
tics of these libraries when N* and iV* are moderately 
large. We now describe these approximations: We as- 
cribe to each state j of the library for a region D(R, N) 
a so-called bulk free energy Q bulk (R, N). Naively this 
would be the sum of the vibrational entropies and in- 
ternal energies of each molecule in D and the pair in- 
teractions between the molecules within D. (Rigor here 
would require ensuring the appropriate finite ranges of 
interactions and some precisely specified ways of parti- 
tioning vibrational entropies: entanglement entropies of 
the interface are too subtle for consideration at present). 
Note that in a system with a globally correlated land- 
scape, such as the harmonic lattice considered earlier, 
a description in terms of local libraries will not provide 
the full list of configurations available to the system as a 
whole. 

In a similar way we would define a bulk free energy of 
the initial state relevant to this region <I>^ fc (R, N). The 
actual energy of the complete sample when the state j 
is inserted in region D(R,N) will not just reflect the 
difference in the bulk energies corresponding to the state 
j and the initial state; instead, it will be higher by an 
amount Tj^n- 

0f - C = (R, N) - $% lk (R, N) + T jtin . (1) 

As mentioned earlier, the new local structure labelled 
by <j> l j h , is likely to be higher in energy than the initial 
configuration. Our construction thus suggests Tj^ n will 
usually be positive and will at most scale with the inter- 
face area. We will write Tj^ n = r yj^ n N x . As we shall 
see later the estimate for the exponent x = 2/3 in 3 di- 
mensions, based solely on the interface area, is probably 
naive near T g . 

We now wish to calculate the equilibration rate of the 
region D(R, N*). If the environment remains frozen, this 
rate is the escape rate from N = to N — N* . The prob- 
ability flux to increase N falls until the bottleneck value 
N* is reached. After N*, even though the average en- 
ergy of states in the libraries increases, the growth of the 
number of states with increasing library size is sufficient 
so that a rapid path to a thermally equilibrated state at 
size N* can be found. Let the typical downhill micro- 
scopic rearrangment have a rate 7"~ icro . This rate will 
only be weakly activated. The flux to any state at a size 
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N smaller than N* will be 



* = r-^J(d0f/^)e s «(»r ,, )/^e-Wr-O/*BT 

(2) 



^3 



r -i P s c ($^"=)/fc B p -(<t>. q -<t>T:)/k B T 

micro ' 



where </> eg maximizes the integrand and is some con- 
stant of units energy. The quantities 4> l j b and $k ulk are 
related through Eq.Q. This maximization means <j) eq 
will be the internal free energy characteristic of the sys- 
tem at the ambient (i.e. vibrational) temperature T. 
That is, 4> eq is the sum of the energetic and vibrational 
entropic contributions appropriate to equilibrium at T. 
In essence, the equation above is nothing more than the 
rate of escape through a transition state with non-zero 
entropy, hence a structure resembling the canonical en- 
semble. Aside from a numerical factor, the integration in 
Eq.pjl is indeed a canonical summation in terms of 4>j, 
which should be regarded as the non-meanficld analog of 
the free energy of the so called "pure" state. The concept 
of the pure state is well developed in the context of frus- 
trated mean-field spin systems 20] . In the mean-field 
limit, the pure states are separated by infinite barriers 
and thus the vibrations around the metastable free en- 
ergy minima are purely harmonic. The quantity (j> l j b is 
non-meanfield because it can be defined only locally (af- 
ter paying the "price" of the surface energy 1^ v ;„). It is 
the degrees of freedom due to transitions between those 
(non-meanfield) "pure" states of the liquid that give rise 
to the (measurable) configurational entropy. Note also 
that owing to the intense and strongly anharmonic mo- 
tions in the liquid (at this finite T), a local liquid state 
labelled by a par ticular value of cj> is not a single inher- 
ent structure |22| , but rather is a weighted superposition 
of many inherent structures, with more or less harmonic 
vibrations on top. 

Note, (j) eq is a function of T and N. Since ^ b ^ lk is 
the equilibrium bulk free energy at temperature T, one 
may replace S c ($ l £f k ) by its equilibrium value at that 
temperature S C (N,T). Thus we get for the escape rate 
to N the result 



k(N) 



expiS c (N,T) 



jib 



k B T 



(3) 



The location of the bottleneck is determined by the min- 
imum of this expression over N . We can thus define an 
activation free energy 



F*{N) = tf> eq -<ffi-TS c (N,T), 



(4) 



whose maximum defines the bottleneck location. Notice 
this activation barrier depends on the total free energy of 
having any target state at size N and the initial nonequi- 
librated particular state free energy <pf b in which we only 
include a vibrational contribution and no configurational 
entropy. This function is shown in Fig[21 Introducing 
Eq.JIJ into Eq.lg} yields: 



F b f k (N, T) - <f> b n k (N, T) + jN x , (5) 




FIG. 2: The free energy to reconfigure the initial configu- 
ration is shown as a function of the size of the locally rear- 
ranged region. There will be fluctuations in the shape owing 
to the detailed packing found in the initial configuration but 
on the average the profile is given by equation in the text. 
The dashed curve shows the profile for an initial state which 
is much higher than the equilibrium energy at T, while the 
solid curve is the free energy profile for an initially equili- 
brated state. 



where F b ^ lk 



(N, T) = -TS C (N, T) is the total equi- 
librium free energy that includes both the configurational 
entropy of the region of size N and the internal bulk free 
energy with its vibrational contribution. In this expres- 
sion we have substituted the average mismatch energy co- 
efficient 7, for the state specific values. For dynamically 
connected states F b ^ tk (N : T), ^ lk (N,T), and S C (N, T) 
grow linearly with N . There will also be 5F, an addi- 
tional fluctuation that typically scales as TV 1 / 2 . Thus we 
write 

FHN) = [f eq (T) - <p in {T)]N + 1 N X + SF, (6) 

where f eq (T) = F b ^ lk (N,T)/N is the total bulk free 
energy per particle of the final state at temperature T 
and (f> in (T) = ^ b " lk (N, T)/N is the internal free energy 
per particle of the initial state. In the following, we 
will omit the fluctuation term SF, whenever computing 
the most probable barrier, but will consider it explic- 
itly when estimating the degree of non-exponentiality of 
the relaxation, which is directly related to the barrier 
fluctuations. Note, the expression above is simply what 
would be prescribed by a regular nucleation theory for 
the free energy barrier of conversion from a (usually non- 
equilibrium) initial state to the other (usually equilib- 
rium) state. Clearly, when the initial state is at equilib- 
rium: & b % lk — & bulk , - the "free energy" driving force 
feq(T) — (f>in(T) is equal to — Ts c , i.e. there is still relax- 
ation in the supercooled liquid, as driven by conversion 
between alternative (local) aperiodic packings of the liq- 
uid. In general, the maximum of the typical F* occurs 
at N* such that dF* /dN\ N=N t vanishes giving 



f, 



eg 



XJ 



and a typical (most probable) rate 



2 

T 



fe 
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(1-x) 



(7) 
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We finish this Section by elaborating on why the re- 
laxations occuring with the rate from Eq.@ can occur 
with a zero free energy gap thus leading to a local free 
energy landscape. Consider first a supercooled liquid 
at equilibrium (i.e. not quenched): /, 



eg 



-Ts r 



Clearly the size N* at which F i (N*) = corresponds 
to the same liquid state, therefore a region of size iV* 
can survey all liquid configurations typical of this tem- 
perature. An isoenergetic state exists because the system 
resides in a metastable state much higher than the lowest 
energy crystalline arrangement, as reflected in non-zero 
configurational entropy s c , and the very high density of 
states e Sc ^\ (Note, this notion underlies the existence 
of residual structural degrees of freedom in glasses even 
at cryogenic temperatures Q.) If the liquid is quenched 
and aging, the driving force f eq — 4>in is still negative 
(it will be computed shortly), so both quantities N* and 
N* exist and are comparable to their equilibrium values, 
although N* no longer signifies straightforwardly the de- 
gree of the landscape locality (consistent with the system 
being out of equilibrium in the first place). Regardless, 
relaxations are local during aging too. 



III. RELAXATION IN THE EQUILIBRATED 
SUPERCOOLED LIQUID 

Examining the expressions for the relaxation rate for 
an equilibrium sample makes explicit the connection to 
the earlier form of the entropic droplet idea. In the equili- 
brated sample case, the internal vibrational free energies 
of the initial and final states are the same so the driving 
force to distinct equilibrated configurations comes only 
from the configurational entropy contribution available at 
the equilibrium temperature. People are sometimes con- 
fused how configurational entropy, which is never made 
available in a strict mean field scenario, can drive a tran- 
sition event. This is because of the locality of the land- 
scape in the droplet analysis, which goes beyond mean- 
field thinking. Locally, for a region of space D(R,N*) 
there is a "funnel" of states leading away from the ini- 
tial configuration to other equally equilibrated configura- 
tions. Configurational entropy drives the activated mo- 
tions of the local regions of a glass in the same way the 
entropy arising from the large number of denatured con- 
figurations of a protein drives the unfolding of a folded 
native protein even though the denatured states are in- 
dividually higher in energy [1|. Any individual escape 
path is likely to find a big barrier but this barrier is par- 
tially cancelled by the growth of the number of escape 
routes. The free energies along each path fluctuate thus 
increasing the likelihood of finding a low energy path, 
when a large number of paths is available. As the en- 
tropy per particle, s c (T) gets smaller the driving force to 
re-equilibrate falls so the critically activated region of the 
glass necessary for re-equilibration grows larger and the 
rate falls, owing to the larger activation barrier from the 
mismatch contribution. The precise way this happens 



depends on the mismatch energy and its exponent x. 

Therefore, before proceeding to study the nonequilib- 
rium situation, we digress to discuss aspects of the av- 
erage mismatch energy T which we have approximated 
as ^N x . First it is clear that the form of jN x is only a 
crude approximation to the mismatch energy when N is 
small. In principle, recall, the mismatch energy could be 
explicitly computed by carrying out the local library con- 
struction on a computer. The only problem is obtaining 
initial configurations equilibrated to the appropriate low 
temperatures. Finding such equilibrated configurations 
deep in the landscape currently requires heroic computer 
resources. Therefore the mismatch energy must presently 
be inferred by analytical considerations. 

The simplest free energy functional calculations of the 
mismatch energy gives an energy proportional to the in- 
terface area and therefore gives the mismatch exponent 
x = 2/3. Such simplified calculations also give an ex- 
plicit value of the prefactor of the scaling relation, 70 
at the ideal glass transition temperature, Tk- At this 
temperature, the lost interaction terms of a sharp inter- 
face have to balance the entropy cost of localizing each 
particle to its cage. 

This gives, at T K HQ, 



7o = — 7: — kbTk m 



{a/d L f 



(9) 



where is the mean square fluctuations of particles in 
a given basin |22l | and a is the interparticle spacing. The 
ratio dh/a, is about 0.1 for glassy configurations, just as 
it is in the Lindemann criterion for melting. Calculations 
based on free energy functional show that 7 vanishes as 
Ta is approached from below Q. This is because Ta 
resembles a spinodal. From these estimates, j(T) can 
be obtained 8]. These naive mean field estimates would 
give an activation barrier varying as s~ 2 , as first detailed 
by Kirkpatrick and Wolynes jj] and later discussed by 
Parisi [28j . 

Kirkpatrick, Thirumalai and Wolynes pointed out an 
effect left out in the naive estimate of the mismatch en- 
ergy p|: One must acknowledge there are numerous so- 
lutions of the mean field equations describing minima of 
the free energy functional, these precisely correspond to 
our local energy landscape libraries. In principle some of 
these other configurations can be interpolated between 
the internal target state to which the region is relaxing 
and its fixed environment in order to lower the mismatch 
energy. This interpolation is called "wetting" . 

Such wetting only can be rigorously defined for very 
large N. Wetting is a dynamical process that takes some 
time to develop. We should therefore for greater accu- 
racy ascribe a frequency dependence to the mismatch 
surface tension y(T,u>). The mapping between the free 
energy functional and the random field I sing model al- 
lows us to invoke an argument of Villain |2!| that gives 
a curvature dependence to the surface tension coefficient 
7 — 7o(o/-R) 1 ^ 2 m three dimensions where a is the micro- 
scopic length where the surface tension 70 is established 
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(R d oc N). Accounting for this wetting correction leads 
to the exponent x — 1/2. When this exponent is used in 
Eq.lJHJl, we see that the relaxation rate in the equilibrium 
system has exactly the Adam-Gibbs form 

k = r m ] cro e- A ^, (10) 

where s c is the configurational entropy per particle. If we 
assume that s c vanishes linearly at an ideal glass tran- 
sition temperature Tx, this rate agrees with the Vogel- 
Fulcher law. Distinct from the AG argument, however, 
in the RFOT theory the critical size N* scales as s~ 2 not 
s^ 1 and the dynamical correlation length is much bigger 
than the AG picture implies, but if the simple estimate 
of 70 based on the vibrational free energy cost is used, 
explicit values of the typical barrier height are obtained 
in addition to the experimental scaling of rates, as found 
by the RFOT theory @- The surface tension Jo/T g , be- 
ing a logarithmic function of the vibrational amplitude 
in the glass, depends little on the atomic make up of the 
glass. Therefore if we use the Vogel-Fulcher expression 

k = fc cxp{-DT g /(T - T K )}, (11) 

the D is predicted to depend inversely on Ac p - the 
change in heat capacity upon vitrification. Specifically 
one finds 

D = 32k B /Ac p (12) 

The resulting correlation of barrier heights (measured by 
D) and glass thermodynamics is excellent [6(. The soft- 
ening of 7 near Ta also explains very well the deviations 
from the VTF laws that are observed as the temperature 
is raised to Ta where there is a crossover to collisional 
dynamics We will use the same "wetted" form of 
the surface energy term without softening in the body of 
the paper to follow. We must note however even at T g , 
these microscopic calculations give N* only of the or- 
der of hundreds of particles (consistent with experiment 
|30j | ) , so we are far from asymptopia and other "ultimate" 
scalings (closer to Tk) are conceivable. 

The local library and mismatch energy concepts let us 
discuss en passant various "defect" pictures of glassy dy- 
namics in landscape terms. Many theories of the glassy 
state imagine there is a basic undefected structure at the 
heart of the phenomenon (an "ur-structure" ) . Doubtless 
for many systems, the periodic crystal itself is one such 
basic "ur-structure". For big enough N* the periodic 
crystal will indeed be a member of the local landscape 
library but it will not be entropically favored and as long 
as the surface tension between liquid and crystal remains 
sufficient it will not be a major target state. Certainly 
"devitrification" can and does occur in the laboratory 
but we will leave the study of this transformation for fu- 
ture discussion. Other ur-structures have been discussed 
as dominating glassy dynamics such as icosahedral crys- 
tals [3lJ or other "avoided crystalline phases" [32] • Deep 
in the local energy landscape, these structures and their 



defected forms must be found. The dimension of the de- 
fects supported by these ur-structures will determine the 
mismatch energy and thus how the rate will depend on 
the driving entropy, in the present picture. Point de- 
fects such as interstitial-vacancy pairs have an energy 
cost independent of the region site i.e. giving an expo- 
nent x — and would give a relaxation rate independent 
of s c . NonArrhenius behavior based on point defects usu- 
ally must rely on special kinetic constraints, which must 
be encoded in the transition rules |33| . By construction, 
the quasi-equilibrium estimate of the rate made above 
would fail for such models. Of course if there were only 
a few "dead-end" states these could be explicitly sub- 
tracted out in the estimates. Such a situation may apply 
for entangled polymers. 

Simple estimates of the free energies of such point de- 
fects puts their energy cost near to the limit seen in or- 
dinary laboratory glass transitions |34j . At low tempera- 
ture they may conceivably short-circuit the generic tran- 
sitions proposed here. Line-like defects enter prominently 
into the constellation of approaches based on frustrated 
icosahedral order [3lJ and presumably should occur also 
when other types of frustrated regular ordered systems 
act as ur-structures. Nussinov has argued that hamil- 
tonians for uniformly frustrated systems should exhibit 
random first order transitions in the mean field approxi- 
mation [3!|. It is likely correct to view the states acces- 
sible in the high density of states region relevant for real 
glasses in this way even for models based on frustrated 
order. 

Recall, in assessing the relevance of defect based mod- 
els, that the measured configurational entropy per parti- 
cle at the laboratory T g is about 1 fc^. This means we 
are usually far from the low defect density regime. 

At low defect density, nevertheless, line-like defects 
may either traverse the correlated region directly or wan- 
der across the region like a Brownian path. The line 
length in the first case of "ballistic" traversal should scale 
like TV 1 / 3 while in the Brownian case the length is pro- 
portional to square of the traversed distance so we expect 
the mismatch energy to scale like TV 2 / 3 , resembling the 
mean field result without wetting. It is interesting that 
an interpolation between the ballistic and Brownian scal- 
ings would be hard to distinguish from the result we use 
r ~ A 1 / 2 . Even at T g , N* is only about 150 so all these 
power laws (with the exception of the point defect case) 
could likely be fitted to a detailed mismatch construction 
with similar accuracy for liquids in the laboratory. 



IV. RELAXATION IN THE IMMEDIATELY 
QUENCHED GLASS 

The microcanonical local landscape arguments allow 
us to estimate the relaxation rates once the statistics of 
the energies of the local regions in the initial nonequi- 
librium state are known. Once the system has fallen 
out of equilibrium, in general, these statistics depend 
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on the detailed quenching history of the sample. The 
complete aging theory should determine these statistics 
seli-consistently. In the present and following sections we 
will assume the quench involves simple straightforward 
cooling and that little time has elapsed since the glass 
transition was passed. In this case the statistics of the 
initial energies are taken to be those of an equilibrium 
system at a temperature T g . Presumably T g will be tem- 
perature of the midpoint of the dynamic heat capacity 
drop upon cooling. Thus we take <^ n (T) to be the bulk 
energy at T g augmented by the vibrational part of the 
free energy at the ambient temperature T. The fluctua- 
tions of the internal energy of regions, which determine 
the nonexponentiality of relaxation in the glass will also 
be taken to be the same as the equilibrium fluctuations 
at T g , and therefore to be determined by Ac p (T g ). 

Using the mismatch exponent x — 1/2 in equations Q 
and (jHJ we obtained fairly simple results for the critical 
cluster size 



XJ 



and the typical rate 
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4fc B T[0 g - f eq ] 



(13) 



(14) 



The initial and final states are at the same vibrational 
temperature but not at the same configurational temper- 
ature. The initial state is typical of a T g configuration 
but the activated state is "equilibrated" to the temper- 
ature T both in vibrational terms and in locally config- 
urational terms. Because of this the driving force for 
reconfiguration is not purely entropic in the nonequili- 
brated case Thus there is a big conceptual difference from 
the usual extrapolation of the Adam Gibbs formula. In 
addition, the rates predicted by Ea. l|14f> differ from the 
Adam-Gibbs formula extrapolated with a fixed configu- 
rational entropy s c (T g ). 

At T g , clearly the rates predicted by the equilibrated 
formula l|10|) and Eg. (|14fl are the same. Upon further 
cooling the nonequilibrium rate while decreasing, how- 
ever, is substantially higher than it would be at equilib- 
rium. This is because the driving force for reconfigura- 
tion not only includes the configurational entropy at the 
ambient temperature but also an energy increment of the 
initial state Ae = e g — e eq (T). (We assume vibrational 
free energy contributions are nearly the same in inher- 
ent structures typical of T g and the ambient temperature 
T.) Just such a change of slope is observed in the labo- 
ratory: the apparent activation energy in the nonequilib- 
rium glassy state is smaller than the extrapolated value 
from the equilibrated supercooled liquid. This change is 
usually quantified in the Nayaranaswany-Moynihan-Tool 
framework through the expression 
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where E* is the equilibrated apparent activation energy 
at T g and xnmt lies between and 1. To compare with 
our initial quench result we would take Tf = T g , again 
noting in a realistic cooling history Tf would need to be 
self-consistently determined. We note that equation 114|) 
gives a gradual transition to Arrhenius behavior. We can 
find xnmt most easily at low temperatures (near Tk), 
where the typical relaxation rate will follow an Arrhe- 
nius law according to Ea. l|14|) . The Arrhenius behavior 
applies because, for ordinary liquids near to Tk we ex- 
pect the mismatch free energy to largely be energetic so 
we take it as a constant. By the definition of the Kauz- 
mann temperature feqiTx) = 4>k is also mostly energetic 
since the configurational entropy vanishes. The vibra- 
tional components of f eq and <p g are assumed to cancel 
so the activation energy is in the nonequilibrium very low 
temperature regime 
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where e g is the energy per particle of the frozen glassy 
state, as prepared, and ex is the energy per particle in 
the putative ideal glassy ground state equilibrated at the 
Kauzmann temperature. The energy difference e g — ek 
is determined by the configurational part of the heat ca- 
pacity Acp(T) for intermediate values of T such that 
T g > T > Tk- The apparent activation energy in the 
nonequilibrated glass at low temperatures turns out to 
be comparable to the equilibrium activation free energy, 
A-F* at T g . It therefore should not vary much from sub- 
stance to substance but depends on the quenching time 
scale tg through the relation 
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To obtain this result of the near equality of A-FJ and 

AE^ L T let us take the configurational heat capacity 
to have the form Ac p = Ac p (T g )T g /T, as suggested by 
Angell |36l l37| . This form is based on good laboratory 
estimates. We can now find by integrating and insert this 
into Ea. (|16fl to give: 
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The activation free energy at T g , on the other hand is 
obtained by finding s c from the integration of Ac p /T 
and using this in the equilibrium rate expression to give: 



AF f = 

9 AAc p {T g )k B T g {T g /T K 
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The ratio of the activation energy in the glass to the 
equilibrated free energy barrier at T g is therefore 
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If Tk and T g are close, as they are for "fragile" systems 
expanding the logarithm gives a ratio close to one. Even 
for the strong liquid Si0 2 the ratio of T g (1480K) to T K 
(876K) is only 1.7, which gives AE^ leLT /AF^ = 0.8, 
if we neglect any temperature dependence of 7. For a 
laboratory glass transition on the one hour time scale we 
should find universally to a good approximation (1 hour 
quench) 32 to 39 k B T g . To compare with the NTM 
phenomenology we note that the RFOT theory predicts 
universally that the activation energy (not free energy!!!) 
in the glassy state depends on the quench time 
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with a coefficient a very close to one. This is consonant 
with many experiments. 

The apparent equilibrium activation energy at T g in 
the liquid is much larger than the nonequilibrium value 
found in the glassy state (being in a sense cancelled in the 



rate expression by a large positive entropy of activation) . 
Again assuming Ac p depends on temperature, we find 
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(22) 



Using the condition that the configurational entropy van- 
ishes at Tk, this energy can also be expressed in terms 
of the ratio of c p (T g ) to s c (T g ): 



(23) 



We then obtain for the inverse of the nonlinearity pa- 
rameter xnmt at low temperatures from Eas. (|20|) and 
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Acp(T g ) I s c (T g ) is a thermodynamic measure of the liq- 
uid fragility. This ratio is large for very fragile liquids 
and small for "strong" liquids. Thus we see "fragile" liq- 
uids are very nonlinear while "strong" liquids, in general, 
should not be. Just such a correlation has been discussed 
by McKenna and Angell [3^ and can be expressed as a 
relation between the kinetic fragility parameter 
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and the nonlinearity parameter xnmt- Using the pre- 
viously derived relation l|20|l and the equation above we 
find 
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Assuming the 7 ratio is near one and using the glass tran- 
sition temperature appropriate to 1 hr, so AF^ /ksTg = 
lnlO 17 ~ 39, and a generic T g /Tx — 1-26 Q, we obtain 
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(27) 



This relation is plotted in Figure [21 along with data for 
several systems. We see that the estimate agrees rea- 
sonably well with experiment. In this estimate we have 



neglected the T dependence of 7. The microscopic treat- 
ment from RFOT shows however that 70 depends both 
on the proximity to Tk and to Ta, as discussed in our 
paper on the barrier oftening effect ||. Thus we will 
generally have both an entropic and an energetic contri- 
bution to 7 o which may explain some of the scatter in the 
curves. Indeed we see more fragile systems lie systemati- 
cally above the curve as is expected since Tk and Ta are 
closer leading to the longer barrier softening effect. 

According to the RFOT theory the typical relaxation 
time in the nonequilibrium quenched state does not im- 
mediately become Arrhenius in temperature dependence 
below the laboratory glass transition. Thus there should 
be deviations from the NMT formalism, which might 
be crudely fit by allowing the nonlinearity parameter to 
be T dependent. Indeed some of the scatter in Fig|3 
probably arises also from this cause. Below T g , the 
expected nonArrhenius behavior from RFOT theory is 
much weaker than the divergently nonArrhenius behav- 
ior found above T g . Within the Adam-Gibbs nonequilib- 
rium extrapolation advocated by Hodge and Scherer 
[l6|. the RFOT theory result would appear to involve a 
temperature dependent "configurational entropy" . Just 
such a behavior has been found by Alegria et al. in poly- 
mer systems |l9l . One way to express this connection is 
to compare AFj e (T) to its value at T g . The inverse of 
this ratio would be the apparent configuration entropy 
in a Hodge and Scherer-Adam-Gibbs extrapolation. In 
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FIG. 3: The fragility parameter m is plotted as a function 
of the NTM nonlinearity parameter xnmt- The curve is pre- 
dicted by the RFOT theory when the temperature variation 
of 70 is neglected. The data are taken from Ref. |3sll - A few 
substances (PVAc = polyvinyl acetate, PVC = polyvinylchlo- 
ride, PS = polystyrene, B2O3, and As2Se3) are labeled. No- 
tice some measured values are not consistent on multiple mea- 
surements; this may reflect a breakdown of phenomenology for 
the history dependence discussed in the text or different ma- 
terial preparation. The more fragile substances lie above the 
prediction without barrier softening, which has no adjustable 
parameters. 



view of Eq. (|14fl , the inverse ratio is given by: 
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(28) 



where the second equality is found by computing f eq with 
the help of the AngelPs form of Ac p (T) (see the Ap- 
pendix). We plot the ratio above in Fig0J We point 
out that although we have highlighted the connection 
with more traditional treatments of aging phenomenol- 
ogy, neither of the earlier approaches is exactly commen- 
surate with our theory. On the other hand we have shown 
the expression for the nonequilibrium rate is very explicit 
once the average energy of the sample is known. It there- 
fore would not be terribly difficult to use the full expres- 
sion i|14|) in a dynamic treatment with the nonequilibrium 
energy as the indicator of the "fictive" temperature. We 
do not carry out this analysis here because it involves de- 
tailed numerical fitting for each system and the quench 
history of each particular experiment. 
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FIG. 4: We plot the predicted variation of the activation free 
energy versus inverse temperature below the glass transition, 
as the ratio (AF*. e . (T)/AF 9 t )" 1 from Eq.l^SJ versus T/T g . 
Only below Tk will the ratio be strictly constant, implying 
strictly Arrhenius temperature dependence of the relaxation 
rate. The temperature dependence of 7 is neglected. The 
figure is plotted for a fragile material with Tk /T g (1 hr.) = 
0.9, a less fragile one with Tk /T g = 0.7 and a strong substance 
with T K /T g = 0.5. 



V. NONEXPONENTIALITY OF RELAXATION 
ABOVE AND BELOW T g 

Within the RFOT theory the mosaic structure of the 
liquid gives rise to dynamical heterogeneity and nonexpo- 
nential relaxation. The driving force for re-equilibration 
varies from mosaic cell to cell. This leads to a range 
of activation barriers, SF$. In the canonical ensemble 
formulation, using the usual Landau formula, these fluc- 
tuations depend on Ac p (T). Computing these fluctua- 
tions allowed Xia and Wolynes to predict /3 for a range 
of substances 0. The same result can be obtained in 
the microcanonical formulation, however here the origin 
is even more transparent: there is a range of energies of 
the initial configurations for each mosaic cell. Since the 
energy fluctuations also scale with Ac p (T), above T g one 
obtains the same results as Xia and Wolynes previously 
derived. 

We first discuss how [3 varies with temperature below 
T g in a somewhat simplified approximation that makes 
the issues clear. We will assume the energy fluctuations 
are small and the resulting barrier distribution is Gaus- 
sian. We will see later that the results of this analysis 
bound the magnitude of the change expected in /3, even 
when a more accurate barrier distribution is used. 

Let us re- write Ea. (|14fl for AF* now including a fluc- 
tuation for the initial energy per particle 
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(29) 



4[</>m - feq + Hin] ' 

If Sdir, is small, we find the fluctuations in the barrier 
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Since the structure is frozen at T g the typical fluctuation 
in (f> is the same as at T g and is determined by Ac p at 
T g . We see the ratio of the size of the AF$ fluctuations 
in the frozen, cooled state 8AF$ e to those found at T g , 
6 AF$ is 
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Both factors in this expression increase rather slowly be- 
low T g , and saturate at Tk- 

How do these fluctuations translate into stretching ex- 
ponents? Roughly speaking, the relaxation function for 
a Gaussian distribution of barriers is approximated by a 
stretched exponential with a value given by 0, E3 : 



0~ [l + (8AF*/k B T)] 
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(32) 



If the liquid is "strong" , Ac p is small so there are small 
energy fluctuations leading to small SAF$. Thus for 
strongly liquids [3 remains near 1 until rather low tem- 
peratures. If the fluctuations in (f> are large (as they are 
for very fragile systems) we would find instead 
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This formula should thus give an overestimate for the 
variation of (3 with temperature. Using this estimate 
along with Eas. (t2$|) and (|3T)l . we find 
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A plot of this expression is shown in Fig OD (note 
the AFj e /AF* ratio has already been computed in 
Eq.||2SJl). At Tk and below the term in parenthesis sat- 
urates and so, we find a simple expression 
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Thus we see in general (3 should fall as we cool but the 
effect remains modest in the range where T is greater 
than Tk- This modestly cooled range is where the most 
detailed aging studies have been reported 0, |4l| . Al- 
though (3 falls, indeed, the rate of fall slows from its 
T dependence above T g for fragile substances, since ac- 
cording to the simple RFOT approximation of Xia and 
Wolynes the equilibrium (3 would vanish at Tk, while 
this (approximate) nonequilibrium vanishes only at ab- 
solute zero. In contrast, (3 for strong substances, that was 
near unity and T- independent above T g , is predicted to 
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FIG. 5: The variation with temperature of the nonequilib- 
rium p in comparison with (3(T g ) is shown as a dashed line 
for a substance with Tk /T g (1 hr.) = 0.8, characteristic of 
a fragile system and a solid line for a strong system with 
Tk /T g = .5. The two other lines indicates how the equilibrium 
would vary with T. The approximte estimates from Eg. 13311 
which exaggerate the variation are plotted as the very thin 
lines, which nearly coincide with the more accurate exression 
Eg. 1321 . The inset shows a magnified view of the region near 
T g , where only factual f3's, depicted by thick lines in the main 
graph, are given. 



show more pronounced tempearature dependence below 
the glass transition (see the inset of FigJSJ. Thus we see 
the form of the relaxation is not the same as the equilib- 
rium relaxation at any "Active" temperature. For fragile 
liquids Tk/T 9 ~ .8, so if 7 is temperature independent, 
a 40% reduction is expected until Tk is reached. Often 
relaxation data in the aging regime have been fit with 
the approximation (3 = (3{T g ). Alegria et al. suggest f3 
remains constant below T g . The RFOT indicates this is 
a reasonable zeroth order approximation. Alegria et al. 
have measured (3 in the regime T g > T > Tk- These are 
difficult measurements and there is scatter in the data. 

We must bear in mind that the theoretical (3 assumes 
that measurements can scan over the complete relaxation 
process (even below T g \) but part of the relaxation is 
missed in experiment. Since the glass only partly relaxes, 
in most experiments (3 will appear to be closer to one 
than if the full relaxation could be followed. This effect of 
missing part of the distribution gives a positive increment 
to (3 at the glass transition, as pointed out by Alegria et 
al. Such an increment on cooling is found. 

After this artifactual increment in f3, Alegria et al. ac- 
tually do find (3 to slightly decrease upon cooling. We 
must emphasize that the arguments leading to Ea. l34|) 
generally give an overestimate of the variation of with 
temperature in the glass state. First and most sim- 
ply the approximation made in going from Eq.(|23) to 
Eq. (l33|) causes an overestimate. Second, but more im- 
portant, as discussed by Xia and Wolynes Q the zeroth 
order Gaussian approximation for the barrier distribution 
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is not quantitative because the simple argument leading 
to that result assumes the environment of a reconfigur- 
ing domain is temporally fixed. Clearly if the domains 
surrounding a region that may re-configure have them- 
selves already changed before the reconfiguration event, 
the library construction's premise of having a fixed en- 
vironment to the mosaic cell fails. This change of envi- 
ronment effect might be called a "facilitation" |3^|. In 
any case this effect means the barrier height distribution 
will be cut off on the high barrier side. A simple cutoff 
distribution follows from the idea that domains slower 
than the most probable rate would actually re-configure 
when their environmental neighbors have changed; thus 
they actually will re-configure at nearly the most proba- 
ble rate, which has already been predicted by the RFOT 
theory. The resulting cutoff distribution for activation 
barrier: 



P(AF*) 



P f (AFi), for AFi < AF * 
C6(AF* - AF|) 



(36) 



- has been shown to reproduce the variation of the equi- 
librium at T g with composition quite well Q . It also re- 
produces the temperature dependence of (3 in the equili- 
brated supercooled liquid. The weight C ensures the nor- 
malization of the distribution. Clearly the cutoff again 
acts to dampen the variation of with temperature, below 
T g . Explicit calculations using Eas. (|36|) and (|31|l for the 
fluctuation can be used to study the detailed T variation 
of (3 in the nonequilibrium glass, which may be relevant 
for deep quenches. 



VI. RELAXATION IN QUENCHES BELOW T K 



bring to mind universal characteristics and indeed the 
RFOT theory does for example suggests a scale for the 
maximum rate corresponding to the cost of overturning 
a single molecular unit; essentially this is 70. On the 
other hand most glass formers have internal structure and 
these inner parts or side-chains of the molecule will have 
multiple conformations that can relax, in a specific ways 
as well those predictable from RFOT alone [42]. These 
"rapid" relaxations, also will slow with temperature, and 
might come into the measurement window. Dyre has 
argued that just such a "contamination" of the main a 
aging process by a (3 relaxation |43j may explain some of 
the Leheny-Nagel measurements on glycerol j4l| . 

2) Well below Tk the mismatch energies may change 
their scaling with size. The "wetting" mechanism re- 
lies on there being a multiplicity of states to interpolate 
through the interface. Such a multiplicity exists at Tk 
but becomes less important with decreasing T. This leads 
to a "hardening" of the interface (increase in 7) or per- 
haps a crossover from iV 1 / 2 scaling to N 2 / 3, scaling of the 
mismatch energy. At the same time, rare fluctuations of 
the barriers may allow some regions to relax more rapidly 
than expected. Plotkin and Wolynes have analyzed just 
this sort of effect in the context of "buffing" of energy 
landscapes in protein Q and Lubchenko and Wolynes 
have studied a similar effect in the quantum regime which 
allows a tail of resonant tunneling states to appear |9(. 
Also well below Tk - small regions may reconfigure by 
crystallizing - an effect we ignore in the present paper. 

3) Eventually below Tk classical barrier crossing will 
be supplanted by quantum tunneling. This gives rise 
eventually to two level systems |9( and the Boson peak 
ll(j which we have discussed in detail elsewhere. 



We have emphasized dynamics in the glassy regime 
just below T g and ranging down to Tk- In this regime dy- 
namics is fast enough so that significant relaxation is still 
accessible to detailed experiments. The results we have 
obtained should hold to considerably lower temperatures 
at least before much aging has actually occured. Essen- 
tially the average rate will become Arrhenius below Tk 
while the breadth of relaxation times will continue to in- 
crease as T approaches the absolute zero. Quantitatively 
several effects may intervene that are worthy of further 
study, however. These effects are 1) secondary relax- 
ations, 2) fluctuations in mismatch energies and changes 
in the "wetting" mechanism, and 3) quantum effects. We 
comment on these in turn. 

1) As a practical matter, less and less of the range of 
relaxation times can be accessed on deep cooling during 
a typical laboratory experiment - at the lowest tempera- 
tures only the fast end of the relaxation time distribution 
can be accessed. This part of the distribution is caused 
by regions of high energy that correspond to a small do- 
main size. At the shorter length scale several effects that 
are system specific, will occur. 

This region of the relaxation time spectrum is often 
called the /3-relaxation. The term /3-relaxation tends to 



VII. AGING AND HISTORY DEPENDENCE 

The local energy landscape theory put forward in this 
paper predicts the relaxation of a system where the statis- 
tics of the energies in the local energy landscape libraries 
of the sample is assumed to be known at any instant. The 
explicit formula for the typical relaxation rate depends 
on the mean bulk energy of the cooperative regions and 
the explicit formula for the stretching exponent /3 also 
contains the variance of these energies. In the quenched 
sample we have assumed these statistics are characteristic 
of an equilibrium system at the temperature T g where the 
system "fell out of equilibrium". Since the cooperative 
regions are large it is natural to assume these statistics 
are Gaussian as we have done. In simple quench histo- 
ries Tg can be estimated by the temperature where the 
apparent heat capacity most rapidly falls during cooling. 
But more complicated thermal histories are possible and 
even in a simple quench T g really must be determined 
self-consistently by the dynamics of the system, itself. 
In general the statistics of the local landscape libraries 
for a nonequilibrium system will be determined by the 
system's dynamics and its detailed past thermal history. 

The most general description of the nonequilibrium 
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FIG. 6: After a considerable period of aging well below T g 
a patchwork of equilibrated and nonequilibrated mosaic cells 
will be found. If the equilibrium energy at T is more than 
a standard deviation of the configurational energies at T g , 
the distribution of energies will be noticeably bimodal and 
the idea of a single Active temperature will break down. The 
unimodal distribution with a single Active temperature should 
be quite safe if AT = T g - T < ^k B T 2 /Ac p N* = ST*. For 
T g relevant to 1 hr. quenches this gives ST* /T g ~ 0.07. Most 
of the Alegria et al. [Tg| data lie in this modest quenching 
range, while "hyperquenched" samples (with AT 3> ST*) will 
often fall outside the allowed range of using a single Active 
temperature. When a sample has a two peaked distribution 
of local energies, an ultra-slow component of relaxation will 
arise. Notice that an equilibrated region at the temperature 
T — T g — ST* will relax on the tens to hundreds of hours scale 
(using the relation that AE^ L T ~ AF%), if r g is taken to 
be one hour. Percolation of the ultra slow regions will lead to 
the possibility of observing multiple length scales. 

statistics is quite complex since the bulk energies of any 
region can be considered functions over the shape and 
size of domains. These functions might be described as 
a set of fluctuating bulk energy fields, but the resulting 
construction is complex. When the system is cooling, 



high energy regions of size N* will be replaced by re- 
gions equilibrated to the ambient temperature T. Thus 
there will be a patchwork of equilibrated and nonequi- 
librated mosaic cells (see Figure |SJ). The nature of the 
new statistics, after some transitions occur and substan- 
tial aging has progressed, depends on how big is the dif- 
ference between the typical current energy of a region 
and the target equilibrium energy. If the gap is big, a 
two peaked structure in the distribution of local bulk 
energies will develop and the statistics will be far from 
Gaussian: some regions that are newly equilibrated will 
relax further at a (slower) rate characteristic of equilib- 
rium at T, while the other, not yet transformed regions 
will still relax at the faster nonequilibrium rate discussed 
already. Such a situation, if it arises, might account for 
"ultra-slow" relaxations which have occasionally been re- 
ported in aging studies H^] (see the Figure caption for 
a criterion for such behavior). The statistics in this case 
of very deep quenches studied for very long times is very 
complex kinematically. Two peaked distributions of lo- 
cal energies may also arise if the system is abruptly but 
briefly heated from a low temperature state to a much 
higher one. Fortunately in the more usual situation of 
modest monotonic quenches we can expect the distribu- 
tion of local energies to remain unimodal. In such cases, 
to first order then the distribution will be characterized 
by a mean bulk energy per particle e! If the ambient tem- 
perature were to remain fixed, e should relax to e eq (T) 
since at that point detailed balance applied to the mi- 
croscopic rates will replenish any states from which the 
system locally escapes. (We remind the reader that e 
denotes a (library) free energy (f> without the vibrational 
entropy contribution.) The typical escape rate at any 
value of e will be given by Ea. (|14fl at ambient tempera- 
ture T. Thus the equation of motion for will e have the 
form 



e-e eg (T) = - j dt^i-t^T^.rteCtO.T^J^tO-e^lTfti)]}. (37) 

J — oo 

I 



where p is the relaxation function when the statistics 
of libraries are known and fixed by e(ir), T(ti) is the 
ambient temperature at time t\ - we assume vibrational 
energies equilibrate quickly by thermal conduction, p not 
only depends on T(t\) and e but contains the parameter 
(3 which in turn depends on the variance of energies of the 
local region. As a first approximation the variance can 
be taken as that characteristic of an equilibrium system 
at a temperature T* such that 



e eq (T*) = e. (38) 



In this approximation, we may use Ea. (|34|) for 

/3 = in [T(t),T(e)]. (39) 

We see that if the shape of the distribution of local li- 
braries does not change much, the RFOT picture leads 
to a situation where a single parameter e suffices to char- 
acterize the system. T*(e) thus essentially fixes the "Ac- 
tive temperature" Tf in the NTM phenomenology, albeit 
with different expressions for t(T,T/) and (3 — (3(T,Tf). 
The non-Gaussian statistics alluded to earlier, however 
suggests the use of a single Active temperature is only 
approximate and that a more complete characterization 
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of the statistics may be needed. Multiple Active temper- 
atures defining the higher moments could in principle be 
defined. In addition spatial fluctuations of Active temper- 
atures are needed to capture the co-existence of equili- 
brated and nonequilibrated domains in the mosaic struc- 
ture. At least for moderate quenches it may be possible 
to ignore the spatial inhomogenity and merely monitor 
the fluctuations in energies of domains as a secondary 
variable. 

VIII. SUMMARY 

We described a local energy landscape theory of the 
dynamics of supercooled liquids and glasses. In the equi- 
librated supercooled regime this theory is just a micro- 
canonical ensemble reformulation of the random first or- 
der transition theory and its notion of entropic droplets. 
New results are obtained in the aging regime of nonequi- 
librium quenched glasses. The key equation is Ijl4(l which 
shows it is the difference between the equilibrated free en- 
ergy at the quench temperature and the initial free energy 
of the particular frozen state that drives motions. The 
theory approximately reproduces the phenomenology of 
Narayanaswamy and Moynihan and Tool. Thus the non- 
linearity parameter in NMT theory can be calculated. 
This parameter is shown to be correlated with the super- 
cooled liquid's fragility in agreement with experiment. 
This correlation is quantitatively very similar to that ob- 
tained by the Hodge-Scherer-Adam-Gibbs extrapolation 
that assumes configurational entropy is fixed at T g . Devi- 
ations from that extrapolation which assumes Arrhcnius 
behavior in the glassy state are predicted however. The 
previously puzzling, modest non- Arrhcnius temperature 
dependence of relaxation observed within the glassy state 
is explained by the RFOT theory, although it is a small 
effect. The variations of nonexponentiality of relaxation 
in the glassy state are predicted but are also rather small 
in the moderately quenched regime. The comparison of 
nonexponentiality with experiment is less conclusive than 
the comparison of mean relaxation rates however, owing 
to the difficulty of accessing the complete relaxation be- 
havior during the quench. 

One advantage of the aging theory based on RFOT 
theory is that in principle the behavior upon very deep 
quenches is predicted. Most importantly (3 is predicted 
to continue to decrease with quench temperature. We 
hope that more experiments in this regime will be done. 
We have noted however that some non-universal effects 
may enter for such quenches. Also the kinematics of deep 
quenches may be complex owing to spatial fluctuations 
of Active temperature predicted by our theory. 

Although the present approach justiAes to some extent 
the use of a single Active temperature to characterize the 
glassy state the limitations of this idea have been made 
clear. A straightforward extension of the usual formu- 
lation to include Auctuations in the moderately cooled 
regime more fully was proposed. In the strongly cooled 
regime, the theory predicts patches of especially slowly 



relaxing regions will appear. This prediction may be 
tested by single molecule imaging approaches j4|| E3] • 
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APPENDIX 

Several auxiliary results are derived in this Appendix 
(here, ks = 1). First, we derive formula (|28|l from the 
main text. Using s = —df/dT, one has (neglecting dif- 
ferences in vibrational entropy) 



feq{T) = 4>k~ 



sdT 



T K 



b K - ACp(Tg)Tg 



T-T K 



K 



where we have used Angell's empirical form |(| for 
the conAgurational entropy, also used in the main 
text: s c (T) = Ac p (T g )T g {l/T K - 1/T), Ac p (T) = 
Ac p (T g ){T g /T). Noting that f eq (T g ) = <f> g - T g s c (T g ) 
Axes the ideal glass state energy 



Ac p (T g )T g ln(T g /T K ). 



(A.2) 



Eas- ljATjl and (|A.2I) immediately yield Eq. l(2"5|) . 

Next, we derive the more accurate expression for the 
j3 n . e .(T)//3(T g ), as follows from Eq.®- We also pro- 
vide the equilibrium value (3 eq of the non-exponentiality 
parameter /3, previously obtained by Xia and Wolynes 
|6(, in terms of experimental parameters Tk/T 9 and 



used throughout this paper. 



First we derive (3 eq . According to Ref.|(|, the bar- 
rier Auctuations are directly related to the local Auctu- 
ations in the conAgurationlal entropy at the scale ./V* 
corresponding to a local equilibrium unit, of which the 
global composite landscape is comprised: SAF+/AF* = 
5S C /S C = y/Ac P N*/s c N*. This yields: 



5AF$ AF* v^cp 



T 



T 



It directly follows from F t (N) = -fy/N 
(dF/dN) N=Nt = 0, and F(N*) = that 

AFt (7/T) 2 



T 



and 



4s, 



7 /T 



(A.3) 



Ts c N, 



(A.4) 



(A.5) 



Using Ea. (|17|l . Ea. (|A.4fl and the temperature indepen- 
dence of the 7/T ratio in equilibrium |6(, one gets 

(7/T) = h/T) T=Tg = 2^s c {T g )\n{T Q /T mlcro ). (A.6) 
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With the help of the equations above and Angell's empiri- 
cal form for the configurational entropy one easily obtains 
that 



Hence, 



SAFf _ y/MTQ/Tmicro) y/l/T K ~ Wg 



T 



T(1/T K - l/T) 



(A.7) 



This and Eq.JJU can be used to compute the 
(3 eq (T)/(3{T g ) ratio. 

The non-equilibrium (3 requires even less effort. Using 
Eci. (|3T|l . one gets 



5AF* 



Fh.(T) 



pi 



l{T) 



(A. 



Pn.eXT) = 



T 



FIAT) 



F 



-1/2 



(A.9) 



The ratio F+ e (T)/F* was computed in the beginning 
of the Appendix and given in the main text as Ea. (|28[) : 
8AF$(T g ) is obtained from Ea.fO|l at T = T g . 
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